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Abstract. A sterile neutrino with mass in the eV range, mixing with v e , is allowed 
and possibly even preferred by cosmology and oscillation experiments. If such eV- 
mass neutrinos exist they provide a much better target for direct detection in beta 
decay experiments than the active neutrinos which are expected to have sub-eV 
masses. Their relatively high mass would allow for an easy separation from the 
primary decay signal in experiments such as KATRIN. 
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1. Introduction 

Currently there is mounting cosmological evidence for an effective number of neutrino 
species slightly larger than three [H El El H]. This result will be tested in detail 
by the Planck mission, and if confirmed will provide solid evidence for physics 
beyond the standard model. Additional sterile neutrino states which are partially or 
fully thermalised via their mixing with the active species is perhaps the most likely 
candidate for extra relativistic degrees of freedom because thermalisation must occur 
well below the QCD phase transition temperature. If not, the density of this extra 
radiation is heavily diluted by the entropy release at the phase transition. 

Intriguingly, also in oscillation experiments there is some evidence pointing to 
the possible existence of sterile neutrinos. The simplest explanation of the LSND 
experiment studying — > v e flavor conversion [5] requires low mass sterile states 
[SI El El El HOI HU H2l H31 HH HSl Uni HZl HH H9] (see also pO] for a more complete list 
of references). 

Very recently, the MiniBooNE experiment has released data from a study of 
the same channel which seems to confirm the LSND result. However, no compatible 
excess was found in the z/ M — > v e channel by MiniBooNE. If both results are correct it 
would require some source of CP -violation. Recent publications suggest two viable 
options; either there are at least two sterile states pTL j |22" 1 [23] . or one sterile state plus 
non-standard charged current -like neutrino interactions (NSI). These interactions 
enables a produced (or detected) neutrino to be a linear combination of flavour 
eigenstates [24J. 

In their paper Akhmedov and Schwetz, [21], find that the 3+1 (active + sterile 
neutrinos) NSI model gives a somewhat nicer fit than the 3+2 scenario mentioned 
above. This is especially true when considering global data - which includes the 
findings of disappearance experiments for v e ^ or v e ^ in addition to the data of the 
already mentioned appearance experiments with — V v e &: — > v e . However 
in the context of beta-decay experiments such as KATRIN [25J (which will be the 
topic of this paper) one measures the number of electrons emitted as a function of 
their energies. So even though the electron neutrino might become a combination of 
flavour eigenstates the electron still has a well-defined energy. Thus, as a consequence 
of simple energy conservation, the electron spectrum is only sensitive to the neutrino 
mass states and a beta-decay experiment would not be sensitive to the NSI. However, 
it should be noted that in the presence of NSI the effective mixing angle appearing 
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in the beta spectrum is not the simply U es , but rather the contribution of the fourth 
mass state to the electron neutrino state, as defined in Eq. 2 of Ref . [24J . 

Considering on the other hand the 3+2 model, the two additional low mass sterile 
states might at first seem incompatible with cosmological neutrino mass constraints. 
But as it turns out the mass bound is actually significantly relaxed by the presence 
of extra energy density, and it is possible to have sterile states with masses close to 
1 eV, compatible with the MiniBooNE result [281 1291 130] . 

Very interestingly in both models the sterile neutrinos would have large mixing 
with v e such that in nuclear beta decay it might be possible to see direct evidence 
for the existence of sterile states. Additionally one would not have to worry about 
CP-phases causing cancellations - as would be the case in neutrinoless double beta 
decay - because the kinematic neutrino mass is an incoherent sum of the mass states. 

In this paper we perform a forecast for KATRIN-like experiments in terms of 
the prospects for observing such sterile states. Even though an active neutrino mass 
of 1 eV is still allowed by beta decay experiments it is disfavoured by cosmological 
data (and would require neutrinos to be Dirac particles in order not to have produced 
observable neutrinoless double beta decays) and we therefore consider a hierarchy 
in which all three active species are very light, with masses too small to be seen in 
experiments like KATRIN, while the sterile state(s) have masses in a range accessible 
by KATRIN (see also e.g. |26l 127] for earlier discussions of this possibility). We then 
investigate how well such an experiment will be able to detect eV mass sterile states, 
depending on their mixing with v e . We wish first and foremost to get an overview of 
the experimental detection probabilities. Therefore we use a rather large sterile 
parameter space in combination with simulations that take all the experimental 
settings of a KATRIN-like experiment into account. The next section contains a 
description of the analysis framework used, Section 3 is a presentation of our main 
results - for both a 3+1 and a 3+2 scenario - and finally in Section 4 we make some 
concluding remarks. 

2. Methodology 

As mentioned we use a realistic simulation of the KATRIN experiment - 
commissioned to begin data-taking in 2012 - to investigate if a similar experiment 
will be able to see direct evidence for sterile neutrino states. 

The main ingredient in the analysis is a toy model Monte Carlo and analysis 
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programme for KATRIN-like experiments written by C. Weinheimer. The 
programme contains all the major experimental settings and has previously been used 
to forecast the experimental sensitivity to the neutrino mass. A much more thourogh 
simulation procedure is currently being built for the KATRIN collaboration. This 
will take into account the finer details of the Windowless Gaseous Tritium Source, 
the differential pumping and cryogenic trapping sections in the transportation of 
the electrons from the source to the spectrometer, and of course the entire final 
electrodynamic properties of the spectrometer. Currently however neither the 
spectrometer or the source have been completed and an ultra-realistic simulation 
might therefore be subject to final changes anyway. So for the purpose of current 
forecasts of the detection abilities of a KATRIN-like experiment this is the most 
realistic simulation tool one can get. For a thorough description of the workings 
and parameters of this code see [37]. We have implemented the x 2 -routine of the 
programme in a modified version of COSMOMC [38j with cosmology turned off and 
flat priors on all input parameters. 

In doing so we have performed a bayesian analysis on the theoretical beta spectra 
while keeping the error bars of the original Monte Carlo driver [39]. It should be 
noted that the mixing angles were held fixed in this analysis. While this is of course 
a perfectly valid way to perform the analysis one may argue that the mixing angles 
should be free parameters as well as the masses. However the mixing angles - or 
more precisely their corresponding entry in the leptonic mixing matrix - determine 
the amplitude of the electron spectra for each separate neutrino mass state. When 
looking at the sum of the mass states, that is the total beta-decay spectrum, it is 
a lot easier to extract knowledge of the neutrino masses (position of the kinks in 
the spectrum) than the mixing angles (relative height of the kinks). So in order 
to get well-constrained results we have kept the mixing angles fixed and one can 
of course perform such an analysis for as fine a grid of mixing angles as desired. 
Finally, oscillation experiments are much better suited for precise measurements of 
mixing angles than beta-decay experiments. In a realistic setting a measurement 
by KATRIN will be combined with short baseline oscillation data to obtain good 
constraints on both mass and mixing angle. 

To the theoretical description of the single neutrino beta spectrum we then add 
massive sterile neutrinos connected to the active neutrino by mixing angles. We 
note that KATRIN can in fact not resolve the mass squared differences between the 
known active states (of Am 2 2 = 8 x 1CT 5 eV 2 and Am^ = 1 2.6 x 1CT 3 | eV 2 ) so in 
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the following we can safely keep treating the active sector as one neutrino^. 

In complete analogy with Eqs. 16 and 21 of [36J the total spectrum is a weighted 
sum of the spectra for the individual mass states. Each mass state is weighted by 
the relevant mixing angles as explained above. That is the appropriate element of 
the corresponding leptonic mixing matrix, U. In the case of the 3+2 scenario the 
full 5-neutrino description U would look like this 
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and a 4 x 4 analogy would play the same role in any 3+1 model. For the remainder of 
this paper the weighting factor is written \U e i\ 2 (V representing the electron neutrino 
and 'i' representing the contribution of mass state i). 

We begin by considering only one sterile neutrino. Under the assumption that 
the three active neutrinos can be treated as one we investigate firstly 2x2 mixing 
schemes where one active and one sterile species are considered as two separate blocks 
connected via a single mixing angle 9. In order to satisfy current mass bounds on the 
active neutrinos we shall assume the active neutrino to be massless and the sterile 
neutrino to have a given mass and mixing angle. We do not investigate the reverse 
mass ordering and we keep the assumption of active massless and sterile massive 
neutrinos for the investigation of the 3+2 scenario as well. In our 2x2 mixing 
matrix we must now have \U ee \ 2 = 1 — |t/ es | 2 , in order to conserve probability. 

We investigate the response of a KATRIN-like experiment to the presence of 
sterile neutrinos for a grid of 11 x 11 combinations of sterile neutrino masses and 
mixing weights. As a starting point for our choice of the parameters m s and U es we 



look at the analysis of 3+2 scenarios in |23] and 
fit area for 2 additional sterile neutrinos: 
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They give the following best 

(2) 
(3) 



f Obviously one could also imagine a scenario where both active and sterile species are massive. 
However we chose to keep the active state(s) massless in this investigation both for simplicity and 
to satisfy current cosmological bounds as far as possible 
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where \U ei \ 2 = sin 2 0(p ->v e ) t i and \U e5 \ 2 = sin 2 Q{p -^v e ),2- Here we want to investigate 
the sensitivity of a KATRIN like experiment to a relatively wide range of mass 
of mixing parameters, including the values above. We take the sterile mass to be 
larger than the projected KATRIN sensitivity of 0.2 eV, and to be conservative 
we investigate masses up to 6.4 eV. We note that such large masses are strongly 
disfavoured by cosmology if the sterile neutrinos are equilibrated in the early universe. 
However, even for large mixing angles, equilibration can be blocked by the presence 
of a e.g. a small but non-zero lepton asymmetry and therefore the cosmological mass 
bound does not necessarily constrain the upper mass range. Meanwhile U es is varied 
from 5.5 x 10~ 4 up to 0.18, i.e. a very wide range of mixing angles. Large values of 
U es are disfavoured by disappearance data (the most recent analysis of which can be 
found in [21]), but are included here for completeness. The exact values of m s and 
U es in the 11 x 11 matrix of parameters studied are listed in Table [TJ 
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Table 1: Sterile neutrino parameters 



We will be calculating the statistical uncertainty on the neutrino masses. 
KATRIN's systematical error is currently estimated to be around 0.017 eV, and 
we shall keep that as our systematic error throughout this investigation. 
In the standard one-active-neutrino case three factors affect the statistical 
uncertainty of the neutrino mass. Firstly the signal count rate at the beta-spectrum 
endpoint. This can be calculated as a combination of the column density of the 
source, the magnetic fields of the spectrometer and source (these determine the 
allowed opening angle), and the cross section of the spectrometer. In our case the 
count rate near the endpoint is controlled via an amplitude factor that is multiplied 
onto the final spectrum (meaning the spectrum after we take final states into account 
and convolute it with the experimental response function). Such a factor - which we 
from now on will call simply the amplitude - depends of course on the experimental 
settings. In the case of KATRIN those settings would be mainly the source column 
density, pd = 5 • 10 17 molecules/cm 2 , the allowed opening angle, 50.77°, and the 
diameter of the effective analysis plane of the spectrometer, 9 m. Given these settings 
the amplitude has the value 477.5 Hz. For reference we include the full calculation 
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of the amplitude parameter in Appendix A. 

The second important issue is the background count rate which can severely 
obscure potentially interesting features near the endpoint of the spectrum. And 
thirdly the energy resolution which determines how well we can see the shape of the 
beta spectrum. 

To sum up: The KATRIN experiment has an amplitude of 477.5 Hz, a 
background rate of of 0.01 Hz and an energy resolution of 0.93 eV. These are the 
parameters we can later tune to simulate an improved Tritium beta decay experiment. 
A full list of the parameter values we have investigated can be found in Table [2j 

For any experiment we want to calculate the detection potential for a second 
neutrino, and the statistical uncertainty on the massless component. 
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Energy resolution [eV] 


0.93 
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Table 2: External parameter settings 



3. Results 

Let us start by examining the sigma detection potential for the sterile neutrino 
component in the standard KATRIN-like setup. Figure [4] in Appendix B shows both 
the detection potential for the sterile neutrino and the standard deviation on the 
mass of the active neutrino. As one would expect, effects are clearest in the sterile 
sector (for which we vary the parameters), but both graphs have a curious feature 
around m s = 1.0 eV. At this point in of parameter space the detection potential 
drops and the standard deviation grows. Inspecting the COSMOMC generated 2D 
likelihood distributions one can see easily that the two mass states are in fact anti- 
correlated up to around ~ 1 eV. An illustration for a mixing weight of 0.18 is shown 
in Figure [U From a simple Taylor expansion of the beta spectrum: 
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we can get the functional shape of the contours in the 2D-plot of Figure [TJ Here 
E e is the electron energy and Eq is the theoretical Q-value (with neutrinos assumed 



to be massless). With the 



m a ^ l_|t/ es |2 Hi's 



-term representing the total neutrino 



mass it is clear that one can write the following expression for the active mass: 

< = m 2 tot - (0.18/0.82) ■ m 2 = m 2 a (m 2 s = 0.0) - (0.18/0.82) • m 2 (4) 

When both masses are very low the mass states act in exactly the same manner in 
the spectrum and the x 2 -function cannot tell them apart. However for high enough 
masses, when the Taylor expansion is no longer valid, the situation changes. Then 
the mass states are not as entangled and the well known correlation between the beta 
spectrum endpoint (which is a free parameter in the analysis - see [37] for further 
explanations) and the neutrino mass comes into play. That is, higher values of E Q 
means higher values of m 2 . Because of the mixing with the sterile mass state, higher 
values of E also mean higher values of m 2 and we get m 2 oc \U es \ 2 • m 2 s . This rather 
weak correlation is the cause of the increase in cr stat (m 2 ) for high m 2 and high \U es \ 2 
that can be seen in Figure HI 

But of course the x 2 -function also depends on amplitude and background. And 
so it turns out that the first point of low sensitivity moves towards lower values of 
m s when we start tuning the experimental configuration. In other words, the signal 
to noise ratio is a decisive factor here. 

While trying to understand what else might affect this situation one should bear 
in mind that the theoretical /3-decay spectrum is folded with the electronic spectrum 
of molecular Tritium and with KATRINs response function (describing the energy 
loss of electrons in the source and the transmission function of the spectrometer). 
Those effects are already included in the code so that we may take them into account. 
Thus the problem should not be the many convoluted functions. However, another 
potentially problematic point of the analysis is the fact that the spectrum is a 
continuous function but the measurements are performed at discrete energy intervals 
around the theoretical Q-value, E Q . Actually the measurement time distribution has 
been optimized for the standard one-neutrino case in which the best sensitivity to 
the neutrino mass is reached when the signal strength equals twice the background 
count rate - see [40j . We have omitted this rather complex optimization for our sterile 
neutrino scenario and that could be a part of the problem. There are certainly hints 
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m 2 [eV 2 ] 

Figure 1: The 2D likelihood distributions for the active (y-axis) and sterile (x-axis) 
neutrino mass-squared for m 2 s = {0.0784, 0.16, 0.3136, 0.64, l^^eV 2 . In all cases 
\U es \ 2 —0.18. The color of the hgures show how many times the COSMOMC Markov 
Chain Monte Carlo (MCMC) has probed a specific area of the parameter space and 
the lines indicate the la and 2a likelihood contours. The mass states are clearly 
anti-correlated for the lowest m 2 -values. Also shown is the line m 2 = m 2 (m 2 = 
0.0) — (0.18/0.82) • m 2 . Clearly the two states become interchangeable and dependent 
on the mass of the other state. This line describes the functional dependence of the 
contours in the 2D plots fairly well until the mass of the sterile neutrino is so large that 
the 2 mass states can be separated. 

that some of the discontinuities in our results may be caused by numerical bad points 
stemming from the measurement time distribution. 

As it is, our results show that for the current settings the sensitivity of a 
KATRIN-like experiment to both neutrino states is lowered unless the mass difference 
is sufficiently large. 

Beyond the signature at 1 eV the graphs shows the basic and expected result that 
a 3cr detection (above the red mesh) of the relic neutrino requires rather high values 
of mixing weight and (or) mass. It also shows that for most of the parameter range 
having a second neutrino in the mix actually improves the statistical uncertainty 
on the mass-squared of the active neutrino. This happens mostly when the extra 
dimension of parameter space due to the small admixture of sterile neutrino absorbs 
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a lot of the uncertainty, resulting in a low m s detection probability. 

If we enhance the amplitude of the experiment one can clearly see the effect 
(Figure [5]). However almost nothing happens for simultaneously low values of \U es \ 2 
and m s as compared to the standard case. Given high values of the sterile parameters 
there is quite some development. But it is certainly not a very smooth development 
as a function of the rising amplitude. This could again hint at numerical discrepancy 
between the 'real' and the measured value of certain parameters as a consequence 
of the measurement time series. In Figure [5] one can also see how the statistical 
uncertainty on the active neutrino gets better for each step - barring the point of low 
sensitivity. This effect definitely corresponds with predictions. By looking closely at 
the figures it can also be made out that the point of low sensitivity moves slightly 
toward lower mass values as the amplitude grows. 

The second important factor is the background count rate. The effect of a lower 
background is shown in Figure [6] of Appendix B and to some extent resembles the 
case for larger amplitude. The development is somewhat smoother and the effect 
not so large, but again there is an enhancement of the detection potential for high 
values of the sterile neutrino components. And as before the lower background has 
an improving effect also on the statistical uncertainty of the active neutrino mass. 

Finally if we try to enhance the energy resolution of the experiment, almost 
nothing happens. Looking at Figure [7] there are only tiny differences in the detection 
potential for the sterile neutrino as compared to the standard case (again only for 
high m 2 s and | C^es | 2 ) ; and apart from that not even the statistical uncertainty on the 
active neutrino mass responds very much. The reason that we get so little effect from 
an improved energy resolution is the molecular Tritium final states. The width of 
the final state distribution is currently comparable to the energy resolution (FWHM 
~ 0.7eV) and so one would have to first avoid that smear of the spectrum (i.e. use 
an atomic source) in order to really gain anything from a better resolution. 

Let us now take a look at an expanded toy model with two massive sterile 
neutrinos to see if another added mass state has any damaging effect on our 
conclusions for the detection abilities of KATRIN-like experiments. Again assuming 
the active states can be described with one single neutrino the 3+2 scenario gives us 
four possible mass orderings (the sterile neutrinos being heavier or lighter than the 
active state, or one light sterile neutrino one massive active neutrino and one more 
massive sterile neutrino). As before the case of two heavy sterile neutrinos and a 
massless active state produces the lowest sum of masses and would therefore be most 
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compatible with cosmological bounds. We did a COSMOMC analysis for this model 

using the best fit point of Eq.'s [2] and [3j The results are shown in Figure [2j 

In the analysis we fixed the values of the mixing angles and again let COSMOMC 

analyse a theoretical KATRIN-like spectrum albeit with Monte Carlo generated 

errorbars. 
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Figure 2: The analysis of a 3-neutrino case using the parameter values of Eq. 's \^and 
The contours (dotted lines) mark the likelihood function in the 2D (ID) -distributions, 
while color (full lines) indicate the data point distribution. The results have converged 
very well on the input values. 
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The output values of our analysis ardjj. 

m 2 sl = 6.47 ± 0.30 eV 2 
m 2 sl = 0.88 ± 0.090 eV 2 
m 2 a = 0.005 ±0.014 eV 2 , 

in agreement with the input values. We observe that our KATRIN toy model 
experiment are certainly able to see also two sterile mass states if such a scenario is 
in fact realized in nature. 

4. Conclusions 

We can see from these investigations that for the highest coupling strengths and 
masses our KATRIN model will definitely be able to separate one or more sterile 
neutrino components from the active neutrino component, if they do in fact have a 
mass and mixing angle in the range proposed by j22], j23], |24| . and [20]. However 
even if just one of the parameters are large there is a very good chance of detection. 

If one wishes to improve the current detection probabilities, the safest way to go 
is by enhancing the amplitude. However it is a rather complex operation to do so. If 
one for instance tries to simply increase the amount of source material in KATRIN 
there will actually not be much to gain. This is because the amount of Tritium 
has been optimized for the size of the source which in turn has been optimized for 
the largest realistic spectrometer-size at the time of construction. Put differently; 
adding more Tritium to the source would not increase the effective column density 
linearly - one would only see a further saturation as illustrated in Figure EJ If one 
then constructed a source with a larger cross section one would also need a larger 
spectrometer. Consequentially - if one used the same strategy as with the current 10 
by 25 m spectrometer - certain smaller german villages would have to be rearranged 
in order to transport such an apparatus to Karlsruhe. 

Another promising option is a lowering of the background count rate. Given 
the numerous background sources in the KATRIN experiment and the continued 
occurence of new ones this may not be feasible either. Obviously the KATRIN 



§ The means are taken from the best fit sample model and the standard deviations from the 
marginalization procedure 
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column density pd [1 15 cm 2 ] 



Figure 3: The ratio of effetive to free column density as a function of the source column 
density in the currently proposed source. Clearly the ratio starts saturating at around 
Pd ~5-10 17 molecules pr cm 2 . The symbols depict the settings for the Troitsk, Mainz 
and KATRIN experiments, and the lines indicate different maximum allowed starting 
angles for the spectrometer. KATRIN has 9 ma x=51 . This figure is taken from the 
KATRIN Design Report J3/J. 

collaboration is making an enormous effort to deal with this problem and maybe 
when the time comes for data-taking it will be easier to assess the opportunities for 
a lowering of the background. 

Finally improved energy resolution will do very little to help these particular 
investigations. However the chances of improving the experimental parameters may 
be better in other future beta decay experiments e.g. MARE or the Project 8 proposal 

For low values of both sterile parameters it looks like not much can be done to 



Detecting sterile neutrinos with KATRIN like experiments 



14 



enable a detection of sterile neutrinos with such properties. At least not within a 
few orders of magnitude of the tunable experimental parameters. 

So in conclusion our results show that given the existence of massive sterile 
neutrinos coupled to the electron anti- neutrino, with the parameters as in Eq.'s |2] 
and [31 KATRIN should definitely be able to see them. In fact for a large enough 
coupling, KATRIN should be able to see a sterile neutrino of any of the mass states 
we have used. Likewise, for a state massive enough, KATRIN would be able to see 
the sterile neutrino for any of the mixing angles we have used. Specifically, in our 
2x2 mixing scheme the experiment could perform a 3<r detection of any of the mass 
states for \U es \ 2 > 0.055. Likewise, for m s > 3.2 eV a 3<j detection could be made for 
any of the mixing angles. In case of la detection \U es \ 2 must be > 0.018 and m s > 
0.8 eV 

Furthermore it is worth repeating that because beta-decay experiments measure 

n 

the incoherent sum: m 2 (z/ e ) = \U 2 J\rn\, any CP phases causing the anti- 

i=0 

neutrino/neutrino asymmetry in MiniBooNE results will not be able to cancel the 
potential signal in beta-decay experiments. Therefore KATRIN is in fact an ideal 
tool for making a direct observation of sterile neutrinos with the properties described 
in this paper. 
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Appendix A: Calculating the 'amplitude' for a KATRIN-like experiment 

The count rate near the endpoint of the Tritium beta spectrum can be written: 
^f(Ee) = N tri ■ e tot ■ K ■ F(E , Z) ■ Pe ■ (E e + m e ) 

n 

Here p e , E e , and m e are the momentum, kinetic energy and mass of the electron and 
m v is the mass of the neutrino. Wi and Vi are the probability and excitation energy 
of the electronic final states. We will now go through the remaining factors of the 
expression one by one: N tr i is the number of Tritium molecules that can be seen by 
the spectrometer: 

Ntri = Pd ■ A S ■ 2 ■ E t ri- 

In this equation pd is the source column density, e tr i is the Tritium fraction in the 
source (that is the source efficiency) and As is the cross section of the source. The 
relation between As and the spectrometer cross section is: 

As — Apinch ■ B max /Bs 

= A^ ■ Ba/ B max ■ B max /Bs 
= A e J f ■ B A /B max ■ 1/ sm 2 (9 max ), 

where A e ^ is the effective cross section of the spectrometer at the position of the 
analysis plane. Bs, B max and Ba are the magnetic fields of the source, pinch and 
analysis plane respectively and 9 max is the maximally allowed opening angle. Keeping 
the dependence on magnetic fields and 9 max in our expressions we now have: 

Ntri = Pd ■ A e / f ■ BA/B rnax ■ 1/ s\n 2 (6 rnax ) ■ 2 • e tri . 

Next up, the total efficiency or acceptance of the experiment, e tot , is a combination 
of the solid opening angle, response function and detector efficiency: 

etot = Afi/4?r • f res (E e , E ) ■ e det 

= 0.5 • (1 - cos(6 max )) ■ f res (E e , E ) ■ e det . 
Here f res is the normalized response function and Edet is the detector efficiency. 
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Finally the factors K and F(E , Z) is a collection of constants and the Fermi 
Function. 

K = G|cos 2 (^)|M| 2 /(27r 3 (V2vr) 7 ) 

= G| cos 2 (5.55 h/{h/27if)/(2n 3 (h/27c) 7 ), 

with Gp the Fermi constant and do the Cabibbo angle. The value of the Tritium 
matrix element, |M| 2 , can be found in e.g. |47| . 



Parameter 


Value 


K 


1.76 • 1(T 17 s" 1 keV~ 5 


F(E ,Z) 


1.19 


Pd 


5 ■ 10 21 m~ 2 


A eff 


(tt • 4.50 2 ) = 64 m 2 


B s 


3.6 T 


B m ax 


6.0 T 


B A 


0.0003 T 


@max 


50.77 ° 


£det 


0.9 




0.95 



Table 3: Selected KATRIN parameters. More details can be found in J37| f 



In total we have 



and 



^(E ) = Pd ■ Af ■ B A /B max ■ 1/ si- 2 ' 



sin 



2 • e tr 



x 0.5 • (1 - cos(6 max )) ■ f res (E e , E ) ■ e det ■ K ■ F(E , Z) 

n 

x p e ■ {E e + m e ) Wi(E -V,- E e ) y/(E - V { - E, 



ej 2 - ml 



i=0 



(E ) = Amplitude • keV 5 (1 - cos(9 max )) ■ f res (E e , E ) ■ p e ■ (E e + m e ) 



x 



Wi(E -V,- E e ) ^(Eo-Vi-E^-ml. 



i=0 



So, in conclusion our definition of the amplitude is 

Amplitude = keV 5 • K ■ F(E , Z) ■ p d ■ A e J f ■ B A /B max ■ 1/ sin 2 



" &det 
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With the parameter values given in Table [3] we get an amplitude of 477.5 Hz 
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lo 9l0 (mf [eV 2 ]) - " log^flUJ 2 ) 

Figure 4: The upper figure shows the sigma detection potential of the massive sterile 
neutrino for the standard KATRIN-like settings, while the lower shows the corresponding 
statistical deviation (in eV 2 ) on the massless neutrino. The x- and y-axis depicts the 
logarithm of the sterile mass squared and the mixing weight. The red mesh illustrates 
the 3a level and the standard one-neutrino statistical uncertainty of around 0.012 eV 2 
respectively (for this analysis method). As one would expect the mass and the mixing 
weight must be rather high in order to get a good detection of the sterile component. 
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Figure 5: The upper figure depicts a comparison of all the amplitude plots for the 
sterile neutrino while the lower figure shows the corresponding comparison for the active 
neutrino component. Included in these graphs is the standard case (being respectively 
the lowest sheet in the upper figure and the highest sheet in the lower hgure). In both 
hgures one can see the low sensitivity point move to slightly lower mass values as the 
amplitude rises. The sensitivity to active mass component gets steadily better as one 
can see from the nice regular sheets in the low ml, |C/ es | 2 area of the lower plot. The 
improvement on the sensitivity to the sterile neutrino component is not a particularly 
smooth function of the amplitude suggesting (again ) numerical bad points in the grid. 
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Figure 6: The upper figure depicts a comparison of all background plots for the 
sterile neutrino. One sees slight improvements on the sensitivity to the sterile neutrino 
component and mostly in the high \U es \ 2 limit. The lower figure shows the comparison 
for the active neutrino component. In both cases one can see the low sensitivity point 
move to slightly lower mass values as the background is lowered. Also we again see the 
stable lowering of the statistical uncertainty on the active neutrino mass component. 
Included in these graphs is also the standard case. 
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Figure 7: The upper figure depicts a comparison of all energy resolution plots for 
the sterile neutrino mass component. In the sterile component there are some slight 
changes where \U es \ 2 is high. In the lower plots one see that the statistcal uncertainty 
on the active neutrino component hardly changes at all. As before the standard case is 
included. 



